Quantum Dissipation due to the Interaction with Chaos 
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We discuss the possibility of having "quantum dissipation" due to the interaction with chaotic 
degrees of freedom. We define the conditions that should be satisfied in order to have a dissipative 
effect similar to the one due to an interaction with a (many body) bath. We also compare with 
the case where the environment is modeled by a random matrix model. In case of interaction 
with "chaos" we observe a regime where the relaxation process is non-universal, and reflects the 
underlaying semiclassical dynamics. As an example we consider a two level system (spin) that 
interacts with a two dimensional anharmonic oscillator. 



The interaction of a system with its environment is a 
central theme in classical and quantum mechanics. The 
main effects that arc associated with this interaction are 
"dissipation" (irreversible loss of energy) and "noise". 
The latter is due to "fluctuations" of the environmen- 
tal degrees of freedom. On short time scales the main 
effect is the "decoherence" due to the noise. On long 
time scales the interplay of dissipation and noise leads to 
a state of thermal equilibrium. 

The common modeling of the system-environment in- 
teraction is provided by a Hamiltonian Htotai = 7~to + 
H(Q, P;x), where H.q is the system Hamiltonian, x is 
a system observable, and TC(Q, P; x) describes the envi- 
ronment including the interaction with the system. The 
simplest (and most popular) modeling of the environment 
is as a large collection of harmonic oscillators. This is 
known as the Caldeira-Leggett approach 0, . Another 
approach is to use random matrix modeling of the envi- 
ronment 0,0, Q However, in this paper we are inter- 
ested in another possibility, where the interaction is with 
(few) chaotic degrees of freedom. 

In what follows "interaction with chaos" means 
that the environment is the quantized version of a few 
degrees of freedom chaotic system. This should be con- 
trasted with "interaction with bath" where the en- 
vironment is modeled as a large collection of quantized 
harmonic oscillators. Let us regard the environment as 
a "black box" (one does not know what is there). The 
questions that we would like to address are: (1) How 
to characterize the bath in a way which does not assume 
a specific model. (2) Is it possible to distinguish "in- 
teraction with chaos" from "interaction with bath" . As 
we explain below the second question is related to the 
notions of "thermodynamic limit" and "universality" . 

Common models for dissipation assume a "thermo- 
dynamic limit" , which means interaction with infinitely 
many degrees of freedom. In this paper we inquire 
whether few chaotic degrees of freedom may have the 
same effect. This is a question of great practical im- 
portance. Future "quantum electronics" may consist of 
several interacting "quantum dots" . One wonders how a 
coherent process in one part of the "circuit" is affected 
by the "noise" which is induced by the quantized chaotic 
motion of an electron in a nearby quantum dot. In other 
words: one would like to know whether it is possible to 



use in the nano-scale reality notions such as "dissipation" 
and "dephasing" , that are traditionally associated with 
having an interaction with many degrees of freedom. 

As a specific example we consider a two level system 
(spin) that interacts with a two-dimensional anharmonic 
oscillator. This would be the well known "Spin-Boson" 
model 0, if the interaction were with a bath of har- 
monic oscillators. The motivation to deal with this model 
is well documented in the cited literature. 

We also compare with the case where the interaction 
is with a random-matrix modeled environment. In the 
"quantum chaos" literature, and in mesoscopic physics, 
random matrix theory (RMT) is regarded as the "ref- 
erence" case. Any deviation from RMT is called "non- 
universality" , and has to do with the underlying semiclas- 
sical dynamics. In this paper we show that the notion of 
(non) universality can be extended into the studies of 
quantum dissipation. 

The basic parameters that characterize any system- 
environment modeling are listed in Table 1. We shall 
define these parameters in a way which is independent 
of the theoretical modeling of the bath. A common as- 
sumption is going to be that the mean level spacing is 
very small. For completeness of presentation, and for the 
purpose of defining what does it mean "very small" , we 
keep A as an explicit free parameter (note ,6|). 



parameter 


significance 


A <x h d 


environment mean level spacing 


Ab oc h 


environment bandwidth 


T 


environment temperature 


d T 


environment heat capacity 


6 


system energy 


A 


system amplitude of motion 


V 


system rate of motion 


r 


strength of coupling 



Table 1: The various parameters that characterize a generic quan- 
tum dissipation problem. (See text for details). 
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As a leading example we consider a two level system. 
For the system Hamiltonian we write Hq — (7if2/2)<Xi 
where is the Bloch frequency, and 01=1,2,3 are the 
Pauli matrices. One can think of this Hamiltonian as 
describing a particle in a double well potential 0, 0- 
Then it is natural to define its position as x = vers, 
where v is a constant. We assume that the interaction 
with the environment is via this "position" coordinate: 
H(Q,P;x) = \{Pl+P$ + Qi+Ql) + (l+x)QjQl This 
environment can be interpreted as a particle moving in a 
2-dimensional anharmonic well (2DW). In the represen- 
tation \v,n), which is determined by H(Q,P;0) and er 3 , 
the Hamiltonian matrix takes the form 



H, 



total 



E + vB hSl/2 
htt/2 E vB 



(1) 



where E = diag{-f?„} is a diagonal matrix that contains 
the energy levels of the environment, and B is a banded 
matrix (sec 9j). The initial state of the total Hamiltonian 
^(t = 0) is assumed to be factorized as ip ® tp, where 
tp is the initial state of the spin, and tp = \no) is the 
initial state of the environment. It is implicit that we 
average over states with E no ~ E corresponding to a 
microcanonical preparation. As for the spin, we would 
like to consider the standard scenario where the initial 
state is a coherent superposition \<p) = (| t) + I !))/v2. 
The reduced probability matrix after time t is 



(2) 



where M — (Mi, M2, M3) is the polarization of the spin. 
It is most convenient to describe the state of the spin 
using M. In particular we define S(t) — M ■ M as a 
measure for the purity of the spin state. 

We turn to formulate the general case. The interac- 
tion of the system with the environment is assumed to 
be of the general form H.i n t — —xJ- where x and T are 
system and environmental obscrvablcs respectively. (In 
the above example x — t>er 3 and T = —Q\Q\). In the 
absence of system-environment coupling we can charac- 
terize the fluctuations of the observable T{t) by a cor- 
relation function C(t). (We are using here Heisenberg 
picture language). Its Fourier transform C E {u) is known 
as the power spectrum of the fluctuations. The observ- 
able T has a matrix representation (n\J-\m) = —B nm 
where we use the basis which is determined by TC. The 
fluctuations are related to the bandprofile of this matrix: 



2 „ r / E m —E n 

nm 1 Z7T0 I U — 



h 



(3) 



n,. . 2irha 2 ( Tiuj \ / huj \ , 
s 2 ™ 6{ " ] + — R \ ~A ) \ At / (4) 

In the first expression there is an implicit microcanonical 
averaging over the states E n ~ E. In the second expres- 
sion er 2 is the average value |B nm | 2 , taken over the near- 
diagonal matrix elements. The lower cutoff function RQ 



depends on the level spacing statistics 0: It is RQ w 1 
for lo > A/h. The mean level spacing A is proportional 
to h d , where d is the the number of environmental de- 
grees of freedom. Any chaotic motion is characterized by 
a finite correlation time t c . Therefore C E (w) has a a cut- 
off frequency ui c = 2tt/t c . This implies (via Eq.©) that 
B nm is a banded matrix with a bandwidth = fko c . 
The envelope function GQ, with G(0) = 1, describes the 
bandprofile. 

For sake of comparison we refer to the Spin-Boson 
model. The distribution of the bath oscillators is de- 
scribed 0, by a spectral function J(uj), leading to 
C E (u>)=2hJ(uj)/(l-e~ muJ ), where (3 is the reciprocal 
temperature of the bath. The "ohmic" assumption of 
having "white noise" (C(t) with short correlation time 
t c ) for high temperatures is imposed by construction, by 
setting J (to) cx luG(lu/lo c ). This corresponds to a strong 
chaos assumption. However, it should be realized that 
even if we "cook" a harmonic bath that has the same 
C e (oj) as that of a chaotic environment, still there is a 
major difference: In spite of having the same bandprofile, 
the B nm matrix of a harmonic bath is very sparse: only 
states that differ by "one photon" excitation of a single 
oscillator are coupled (hence the differences E m —E n are 
the frequencies of the harmonic oscillators). 

In what follows we would like to define the notion of 
temperature (T) without assuming a specific modeling. 
There are three features of the environment that has to 
do with this notion: (i) The growing density of states; 
(ii) The growing fluctuations intensity; (iii) The asym- 
metry of C e (uj) with respect to uo. Let us regard the sys- 
tem as a "thermometer" . The equilibrium is determined 
by the microcanonical temperature T — (&e ln(l/A)) _1 . 
For d=2 environment with constant density of states we 
would get T=oo. [We consider this hypothetical case for 
argumentation purpose. This would be indeed the case 
if the environment were modeled as a billiard. In a later 
paragraph we discuss the model of Eq. for which T is 
finite]. In case of a two level system, having T = 00 im- 
plies an equal probability for the two energy states. Now 
we can take (instead of a two level system) a particle with 
one degree of freedom as a "thermometer" . Since we al- 
ready know that the "temperature" is (say) T = 00, we 
may deduce that the friction coefficient, as determined by 
the fluctuation-dissipation theorem is 77 = vj(2T) = 0, 
where v — C e (uj^0) is the intensity of the fluctuations. 
This conclusion is wrong. In fact the friction coefficient 
is 7] = iAx3^(^/A), as discussed in Q (and references 
therein). Thus we can define an effective temperature, 
which is related to the thermalization process: 



v 

2^ 



^ln(i(7 E (^0) 



(5) 



In generic circumstances the distinction between T eff and 
T is not so dramatic. In fact, some further inspection into 
Eq.© reveals that generically T etl = T/2. The above 
subtlety does not apply to a thermal preparation. As- 
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suming a reciprocal temperature (3, the friction coeffi- 
cient (as obtained by canonical averaging over the above 
cited microcanonical result) is r\ = \(3v. Therefore we 
get T eff =l//3 as expected. 

Having defined T, we can make a quantum mechanical 
distinction between high and low temperature regimes. 
This is related to the asymmetry of C E (w) with respect to 
uj. Some inspection into Eq.@ reveals that a finite tem- 
perature implies that the band profile acquires a factor 
exp(cj/(2T cff )), which is consistent with the Spin-Boson 
modeling (see a previous paragraph). "High tempera- 
ture" means that for the physically relevant frequencies 
Tiu/T ^ 1. A sufficient condition is T ^> A&. In the 
latter case C E (u>) can be treated as a symmetric func- 
tion with respect to w h -w, and therefore it can be 
interpreted as the power spectrum of a classical noise. 

The issue of thermodynamic limit is related to the heat 
capacity of the environment. This is d T = (dT / dE)^ 1 ~ 
d, where d is its number of environmental degrees of free- 
dom. The energy that the system can exchange with the 
environment is denoted by €. If we want to assume a 
stable temperature T, the heat capacity of the environ- 
ment should be large enough, so that energy exchange 
between the system and the environment does not have 
a big effect. This leads to the condition 6 <C d T x T. 
For a generic few degrees of freedom system € ~ c?o x T, 
where do is the number of system degrees of freedom. 
This leads to do <C d T . In case of a two level system the 
condition is much easier. Namely, we have 6 ~ Ml, and 
therefore we get the "easy" condition Ml <C d T x T. 

Assuming the typical circumstances of an oscillating 
system, the time variation of the system observable x(t) 
is characterized by an amplitude |x| ~ A, and by a rate 
of change \x\ ~ V. This specification of the system dy- 
namics is essential in order to define a dimensionless pa- 
rameter that characterizes the system-bath interaction: 

£ = minimum (^A\ (%V) 2/3 ) (6) 

Loosely speaking this parameter indicates how many en- 
vironmental energy levels are mixed non-perturbatively 
due to the interaction with the system. The A depen- 
dence of r is the consequence of the well known theory 
by Wigner: It is the number of levels which are mixed by 
the perturbation in 7i i— > E n 8 nm + xB nm with |a;| ~ A. 
If the perturbation is slow (small x) this A based esti- 
mate becomes non-relevant. For a proper analysis [g] one 
should switch to the adiabatic (x dependent) basis, lead- 
ing toHn E n 5 nm + x(iKB nm /(E n -E m )) with |i| ~ V. 
The V based estimate for T follows from the latter rep- 
resentation. Thus one realizes that for slow rate of x 
variation there is a crossover form a A-determined to a 
1^-determined mixing. 

The rest of this paper is aimed in clarifying the sig- 
nificance of the parameter T. Some results of the sim- 
ulations with the 2DW model Eq.0J are presented in 
Fig. I. The energy of the environment was in the range 
2.8 < E < 3.2 where the classical dynamics is predomi- 
nantly chaotic. The classical correlation time is r c ~ 1. 



The simulations are done with fi = 0.03. This means that 
the bandwidth is At, ~ 0.2. This should be contrasted 
with the mean level spacing A ~ 0.004. The tempera- 
ture in the specified energy range is T ~ 1.3, and the 
heat capacity is d T ~ 2.4. The amplitude of the motion 
is A = v, while the rate of the motion is formally ^ = 00. 
The latter should be understood in the path-integral con- 
text (Eq.Ql), where x(t) makes "jumps" between the two 
sites (x = ±v). Hence T oc (v/Ti) 2 . One can easily verify 
that the Kondo parameter Q is a = (l/16n)T/T. In our 
simulations a -C 1. 

In the lower panels of Fig. I we present the correspond- 
ing results of simulations with a RMT model where the 
induced fluctuations have exactly the same power spec- 
trum as in the 2DW model. The RMT model has been 
obtained by taking the Hamiltonian Eq.Q with a ran- 
domized B. We simply randomized the signs of the off- 
diagonal elements. This procedure destroys all the cor- 
relations between the elements, but does not affect the 
bandprofile (which is implied by Eq.JSJ)). Two observa- 
tions should be made immediately: (i) A few degrees of 
freedom chaos indeed provides a dissipative effect, as in 
the case of a many degrees of freedom bath, (ii) The 
effect of interaction with "chaos" can be distinguished 
from the case of a random- matrix modeled environment. 
The latter claim is based on the observation that in the 
regime v > 0.16 there is a two orders of magnitude differ- 
ence between the corresponding curves. Also the scaling 
of curves with respect to v is "broken" (upper panel): 
The sensitivity to v is much smaller than implied by the 
overcompensating scaling of the time axis. 

As explained in the introduction the deviation from 
RMT is regarded as "non-universality" and its under- 
standing requires a proper definition of the classical 
limit. This should not be mistaken as a synonym 
for " high temperatures" . The confusing, and possibly 
meaningless procedure to define this limit is by taking 
h — > 0. We argue below that the meaningful definition 
of a "semiclassical regime" is by considering the nature 
of the dynamics of the environment. This leads to the 
condition T At,. In order to explain this condition 
we adopt the Feynman- Vernon picture of the dissipation 
process. Within this picture the propagator of the re- 
duced probability matrix is written as a path integral: 

K(v,v'\v ,v' )= F[xa,x b ] e^^-^" (7) 

X A ,X B 

The summation is over pairs of system trajectories (with 
weight factors absorbed into the definition of the inte- 
gration measure). The action A[x] is defined as the 
phase which is accumulated along a given trajectory. 
We note that for a spin the trajectory is piecewise 
constant (x — ±v). The influence functional is de- 
fined as F[xa,xb] — (^\U[xb]~ 1 U[xa}\'<P} where the 
expectation value is taken for the initial preparation 
of the environment (it is typically a mixture of many 
states implying that an appropriate average should be 
taken over ip). The environmental evolution opera- 
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tor in case of the models that we consider is U[x] = 

T exp (-(i/h) f*(E + x(t')B)dt'Y Thus, in a semiclas- 

sical framework the problem of "quantum dissipation" 
reduces to that of analyzing "driven degrees of freedom". 
It has been realized that the driven dynamics be- 
comes non-perturbative if T > At,. Below we discuss the 
implication of this observation. 

The purity S(t) is related to the dephasing fac- 
tor IF^a, xb}\- [For Q—0, the Hamiltonian Eq.QJ 
becomes block diagonal, leading to S(t) = \F[xA, %b]\ 2 , 
with xa=v and xb=—v\. If the fluctuations are re- 
garded as "noise" , one obtains the standard expression 

pr i| --rhr V f * C(t'-t")(x B (t')-x A (t")) 2 dt'dt" 

\t [xa, xb\\ = e ffl2 JoJo 

This expression implies a short time Gaussian decay 
S(t) = exp(— 4:C(0)(vt/h) 2 ), which evolves into a long 
time Gaussian decay S(t) = exp(— (avt/Ti) 2 ) in the 
regime T < A. If the envelope G() of the bandpro- 
file were smooth, then we would expect, in the regime 
A < r < Af,, an intermediate stage of exponential de- 
cay S(t) = exp(-2 7 i) with 7 = 2(v/h) 2 C E (Q) ~ T/H. 
In case of our numerical example the bandprofile has a 
structure 0, and hence C(r) has oscillations that show 
up in the simulations. 

Can we trust the standard expression for ^[x^, xb]\ 
if we have a dynamical environment rather than a noise 
source [with the same C(r)]? It is not difficult to observe 
that the standard expression can be derived from Fermi- 
goldcn-rule (FGR) considerations. In case of Harmonic 
bath the FGR treatment is valid, and this expression, 
with C(r) replaced by its symmetrized version, is ex- 
act. But for interaction with "chaos" we have claimed 
above that non-universality should show up in the non- 
perturbative regime (r > A;,). Our numerics confirm this 
prediction: In the non-perturbative regime (r > At,) we 



find a premature (t < t c ) crossover from the expected 
short time Gaussian decay, to a non-universal behavior 
which is not captured by the standard formula. In the 
2DW case the system has a classical limit, and there- 
fore the nonperturbative decay of S(t) is slowed down 
(compared with RMT) because it is limited by the clas- 
sical dynamics. As observed (note the inset) the decay 
becomes much less sensitive to v. 

The above discussed non-universality can be regarded 
as the manifestation of "semiclassical" correlations be- 
tween the off-diagonal matrix elements of B nm . These 
are not reflected in C e (c<j). In case of the Spin-Boson 
model the B nm matrix is sparse, which implies (in the 
limit of infinite bath) that off-diagonal correlations can 
be neglected. In case of a random matrix modeled en- 
vironment, absence of correlations is guaranteed by con- 
struction for x = 0, which implies (assuming no sparsity) 
lack of "invariance" with respect to x Accordingly 
we have three classes of models that become distinct in 
the non-perturbative regime. 

In conclusion, we have discussed the consequences of 
having four energy scales (A,At,T, T) in any generic 
problem of "quantum dissipation" . The strength of the 
interaction is characterized by T. A universal "quan- 
tum dissipation" behavior requires a separation of en- 
ergy scales A r <C At,. Non-perturbative breakdown 
of this universality due to the underlaying semiclassical 
dynamics is found if T > A;,. 
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FIG.l: The decay of S(t) in case of interaction with chaos (upper panels) and in case of interaction with a "randomized" environment 
that has exactly the same fluctuations (lower panels). The selected values of the coupling parameter arc in the range 1(T 4 < v < 0.3. The 
left panels are for the T < A regime (v < 0.02) and for the T > A;, nonperturbative regime (v > 0.16). The inset is for some of the v > 0.16 
curves without scaling. The arrows point on the largest v value. The right panels are for the A < T < A;, regime (0.02 < v < 0.16). The 
scaling methods of the time axis are implied by the analysis of S(t), which should hold if there is a universal behavior which is determined 
by C(t) alone (see text). 



